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Abstract 

We obtain explicit expressions for the correlation functions of transmission and 
reflection coefficients of coherent electronic waves propagating through a disordered 
quasi-one-dimensional medium with purely elastic diffusive scattering in the pres- 
ence of spin-orbit interactions. We find in the metallic regime both large local 
intensity fluctuations and long-range correlations which ultimately lead to universal 
conductance fluctuations. We show that the main effect of spin-orbit scattering is to 
suppress both local and long-range intensity fluctuations by a universal symmetry 
factor 4. We use a scattering approach based on random transfer matrices. 
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I. Introduction 

There has been considerable effort made in the past few years to understand the 
statistical properties of the intensity pattern created by coherent waves propagating 
through a disordered medium in the multiple elastic scattering regime. The simplest 
known effect is large local intensity fluctuations (the speckle pattern) observed in 
scattering of light [1], which can be interpreted theoretically [2] as being due to 
fluctuations in the sum of a large number of uncorrelated amplitudes associated with 
different Feynman paths arising from the diffusive motion of the scattered wave. 

Another interesting feature of such patterns, which has also been experimentally 
observed [3,4], is the presence of coherent backscattering. A result interpreted as a 
constructive interference between pairs of time-reversed scattering paths. 

It is now understood that the statistical properties of the fluctuating intensity 
patterns may crucially depend on the complex interference of the underlying wave 
fields of the diffusion modes in the bulk of the medium. This dependence leads 
to non-trivial long-range correlations [5,6], which in turn imply fluctuations of the 
total transmission coefficient of order unity. Such fluctuations have been observed in 
phase coherent electronic waves transmitted across a disordered sample [7-10]. The 
theory of this effect [11-15] predicts that the conductance of a disordered mesoscopic 
metallic conductor fluctuates with an amplitude always of order e 2 /h independent 
of the system size and degree of disorder, the phenomenon of universal conductance 
fluctuations. 

An important consequence of this result is that it implies strong correlations 
between diffusion modes and thus reflected and transmited intensities in different 
modes must be statistically correlated. 

All these interesting statistical properties are, in fact, expected to occur in any 
system in which waves propagate coherently in a disordered medium. Electronic 
waves, though, can give rise to special quantum mechanical effects if the impurity 
scattering becomes spin-dependent by adding to the sample strong spin-orbit cou- 
plers like Au. The simplest observable effect is the depletion of backscattering [16], 
which results from coherent destructive interference of spin-dependent wave func- 
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tions associated with pairs of time-reversed scattering events. Theoretically, it has 
been shown [17] that spin-orbit scattering supresses universal conductance fluctua- 
tions by a factor of 4 and gives rise to weak anti-localization, slightly enhancing the 
value of the average conductance. 

There are two different approaches to coherent electronic wave propagation 
in disordered medium. The first one [11-13] (the microscopic approach) is based 
on linear response theory and uses Feynman diagrams to perform averages over 
different impurity configurations. The second one [15,18] (the macroscopic approach) 
is based on Landauer scattering theory [19] and uses random transfer matrices to 
represent the multiple elastic scattering of electronic waves in the sample. These 
two approaches are essentially equivalent for systems with quasi-one-dimensional 
geometry. 

Explicit expressions for the correlation functions of the transmission coefficients 
of waves propagating through a random medium with purely elastic scattering have 
recently been obtained by Feng et al. [6] using the microscopic approach. They 
obtained, in addition to the familiar speckle pattern, long-range correlations resulting 
from the interaction of diffusion modes. A similar result has also been derived by 
Mello et al. [20] using the macroscopic approach. They showed that long-range 
correlations between reflection coefficients are of the same order as that between 
transmission coefficients, a result which had some importance in connection with a 
deeper understanding of universal conductance fluctuations. 

This paper is concerned with correlations of transmission and reflection coef- 
ficients of electronic waves diffusively scattered in a disordered medium both with 
and without strong spin-orbit coupling. To the best of our knowledge there is no 
explicit evaluation of correlation functions of transmission and reflection coefficients 
in the presence of spin-orbit scattering. In this paper we provide such a calculation 
together with results for a correlation function not considered in previous works. In 
presence of spin-dependent scattering a number of new purely quantum mechanical 
interference effects takes place and the above authors results do not apply, even lo- 
cally, where we found the amplitude of intensity fluctuations to change dramatically. 
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We use the macroscopic approach to obtain results for the unitary, orthogonal and 
symplectic ensembles. The unitary ensemble is appropriate to systems with an ap- 
plied magnetic field. The orthogonal ensemble describes systems with time-reversal 
symmetry in the absence of spin-orbit scattering, whilst the symplectic ensemble 
applies to systems with time-reversal symmetry in the presence of strong spin-orbit 
scattering. In Sec. II we formulate the scattering problem and show how the various 
quantities of interest can be written in terms of some parameters that appear in 
the polar decomposition of the transfer matrix M of the disordered sample. We 
assume the system to be quasi-one-dimensional, which greatly simplifies the calcu- 
lations owing to the isotropic nature of the probability distribution P(M) in the 
group manifold of the transfer matrix. In Sec. Ill we obtain novel exact expressions 
for the various correlation functions of reflection and transmission coefficients in 
terms of the number N of transverse propagating modes and certain averages to be 
performed by solving an isotropic diffusion equation for P(M), in which the length 
L of the sample plays the role of an embedding parameter. In the metallic regime 
we find that large local fluctuations and long-range correlations are suppressed by 
coherent interferences generated by spin-orbit scattering. We show that our result 
for the speckle pattern can be understood in terms of a sum of a large number of 
uncorrelated Feynman paths. The importance of long-range correlations to universal 
conductance fluctuations is briefly discussed. We present a summary and conclusions 
in Sec. IV. 

II. The Scattering Approach 

We consider a confined quasi-one-dimensional geometry in which semi-infinite 
perfectly conducting leads are attached to the right and left sides of a disordered 
region of length L and width W. The total number of transverse propagation modes 
is given by N ~ where kF is the Fermi momentum and d is the sample 

dimension. We assume L <C L in , where L in is the inelastic scattering length, thus 
the only important scattering mechanism is elastic and the electron wave function 
is coherent. 
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Since the orthogonal and unitary ensembles have been thoroughly discussed 
in many recent papers [15,18], we present a brief review only of the less familiar 
symplectic case. 

The amplitude of the propagating modes to the left of the sample are related 
to those to the right by the transfer matrix M. Flux conservation and time-reversal 
symmetry in the presence of spin-orbit scattering imply that M can conveniently be 
parametrized as [17] 

v^Xo :•)• ™ 

where u and v are arbitrary N x N quaternion unitary matrices and A is an N x N 
quaternion real, diagonal matrix whose elements Ai, A2, . . . , Ajv define the eigen- 
channels. The operations over quaternion matrices and elements are defined in Ref. 
17. 

The transfer matrix M can also be written in terms of transmission and reflec- 
tion matrices t, r (fluxes incident from the left) and t' , r' (fluxes incident from the 
right). It reads 

M-( (* f ) _1 r '^~ 1 ] (2 2) 

Comparing (2.2) with (2.1) one finds 

1 

1 \ 2 



t = n {—x) v - (2 ' 4) 

The reflection and transmission coefficients R a b and T ab of wave amplitudes in chan- 
nel a when there is an incident flux in channel b are given by [21] 

iUHrJ 2 = $>2| 2 , (2.5) 

i=0 

r a6 HU| 2 = X>2l 2 , (2-6) 

i=0 
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in which r ab and t ab are complex components of the quaternions r a b and t a b in the 
standard basis {e^, i = 0, 1, 2, 3}, where eo is the 2x2 unit matrix and e±, e% are 
defined as 

We define three kinds of crossed second moment 



C a b,a>b> = (TabTa'b'} ~ (T a b) (T a > b') , (2.8) 

D a b,a'b> = (RabRa'b') ~ (Rab) (Ra' b') , (2.9) 
E a b,a'b' = (RabTa'b') - (Rab) (T a 'b') , ( 2 -10) 

in which (...) stands for an average over an ensemble of random matrices M of 
the form (2.1). It has been shown [17] that, as the sample length is increased, M 
performs an isotropic random walk in its group manifold. This implies that ensemble 
averages of arbitrary functions of the form f(M) = f 1 (u 7 v)f 2 (X) are simply given 
by 

(/(M))H/iM>o</2(A)> s , (2.11) 

where (/i) denotes the average on the quaternion unitary group and (f 2 ) s is the 
expectation value obtained by solving the evolution equation 



o(h), L/lfl 

-l\J^dX a 



ds 2N .. , 

a=l 



A a (l + A a )j|^ 



(2.12) 



in which J = Yl a <b l^o — | 4 and s is the length of the sample in units of the mean 
free path. 

The quasi-one-dimensional geometry owes its mathematical simplicity to the 
factorization shown in (2.11). 



III. The Correlation Functions 

If we insert (2.6) into (2.8), use (2.11) and perform the averages on the quater- 
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nion unitary group in the way described in Refs. 17 and using the results of Ref. 
22, we find 



C ah , a , b , =(A C , 1 (T 2 ) S - A Cj2 (T 2 ) s )(4 + 8 aa ,8 w ) 

+ (Ac,i(T2) s -4Ac,2(T 2 ) s )(Saa'+6 bb/ ) - (T ab )(T a , b ,), (3.1) 

where 

4N 2 + 1 2 
Ml ~ AT2 (4 iV2-l)2' Ac > 2 ~ iV(4iV2-l)2' 

JV TV 
a6 a=l a=l 

Equation (3.1) is exact. As a check one can verify that 

Var(T) = C ab,a'b> (3.4) 

aba'b' 

holds for arbitrary N. 

Inserting (2.5) into (2.9) and performing the averages on the quaternion unitary 
group we get 

D ab , a , b/ =A DA ((R 2 ) S - (R 2 ) s )5 aba , b , + (A Di2 (R 2 ) s - A D , 3 {R 2 ) S )((2 - 5 ab )(2 - <W) 
+ <W<W + 5 ab iS b a>) + 2(A D ^(R 2 ) S - A DA {R 2 ) s )(S aa 'b' + 5 a w + S aba ' 
+ 5 ba > b > — 5 aa > — 5 ab > — 5 bb > — 5 ba >) — (R ab ){Ra'b'), (3-5) 

in which 



A 



D,2 



A D a = 





1 


(3.6a) 


N(N - 


- 1)(2AT-3)(2AT- 1) ' 




2N 2 -3N + 2 


(3.66) 


N 2 (N- 


-1)(2N -3)(4iV 2 -1)' 




1 


(3.6c) 


N 2 (N 


-l)(2iV-3)(2iV+l)' 




2N 2 - N + 1 


(3.6c?) 


2N 2 (N 


- l)(2iV-3)(4iV 2 - 1)' 



and 



N , N , . x 2 



«=E^ = Ei^. * = E(i^ • < 3 - 7 ' 

ah a=l a=l 
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The unit tensors 5 abc and S a bcd take the value 1 if all indices coincide and vanish 
otherwise. One can easily check that D a b,a'b' satisfies 



Var(#) = D *b,a>v. (3.* 



aba'b' 

Finally for (2.10) one finds 

E ab , a >b> ={AeART) s - A Ei2 ((T) s - (T 2 ) s ))(2 - 5 ab ) + (A Efi ((T) a - (T 2 )J 

- A E , 2 (RT) s )(6 ab > + 6 bb > - <W) - (Rab)(T a >b>), (3.9) 

where 

Ae ' 1 = N*{2N+l){N-iy NAE > 2 = Ae * = N(N - 1K4AT2 - 1) ' (3 - 10) 
The averages (R a b) and {T a b) have been calculated in Ref. 17 and are given by 

Observe that (R aa ) < (Rab) in agreement with weak- localization theory [16]. One 
can verify that E abja > b / satisfies 

Cov(R,T)= E ab,a>b>. (3.12) 

aba'b' 

Equations (3.1), (3.5) and (3.9) are new predictions of the present macroscopic 
model. Comparing equations (3.1) and (3.5) with the ones obtained in Refs. 15 
and 20 for the unitary and orthogonal ensembles we find great similarities in their 
structure. We remark that the correlation function E ab , a 'b' has not been considered 
in earlier works. 

To proceed further we have to evaluate the remaining averages using the evolu- 
tion equation (2.12). Since no exact solution is possible for arbitrary s, we confine 
ourselves to the important metallic or diffusive regime, defined as the one in which 
1 C s < f. In this regime the average total transmission (T) s , which is propor- 
tional to the conductance, has a very large value and it makes sense to expand an 



arbitrary average (/(A)) s in inverse powers of (T) 8 . This procedure is explained 
thoroughly in Ref 15. We merely quote here the resulting expressions for C abia > b >, 
D abja > b > and E ab y b >, which are 

(1 2 2 \ 

-Saa'Sbb' + » /rrx (<W + &bb') + .rn/^o ' (3.13) 
a 3a(T) s 15(3{T)jJ 

n /D \/D \( P (l5aa>&bV +l&ab>ha> 

JJ a b,a'b' — \tt a b)\Ka'b' 



2a V /? + <U(2 - (3) 



(2(5 aa > + <W) + l{5ab> + ha' ~ f3Saa>bb>)) + 32 . ^ ^ 



2a(R) s 15[3(R) 



E ab , a , b , - -(R ab )(T a , b ,) {^^f + m R )g{T)s ) ' ( 3 - 15 ) 
where a = (f3 2 - 3/3 + 4)/2, 7 = (2 — /9)(3 — /?) and /3 takes on the value 4 in the 
symplectic ensemble. Other possible values are (3 = 1 (orthogonal ensemble) and 
(3 = 2 (unitary ensemble), which have been studied in Refs. 15 and 20. Expression 
(3.15) is, however, for all values of /?, a new result of the present work. 

Suppose a single mode b is excited on the left and we are interested in corre- 
lations between transmission coefficients in different modes a and a' on the right. 
From (3.13) we get 

Cov(T ab , T a > b ) = -(T ab )(T alb )6 aa >, (3.16) 
a 

thus the present model predicts large local fluctuations of transmission intensities. 
For conventional waves, a = 1, and (3.16) is the familiar speckle pattern observed in 
many experiments [1]. For electronic waves in the presence of spin-orbit scattering, 
a = (3 = 4, and from (3.16) we see that local fluctuations are suppressed by a factor 
of 4. 

The above result can be easily understood by means of a very simple argument. 
The presence of disorder implies that the electronic motion is diffusive and thus the 
transmission amplitude t ab is a large sum of amplitudes associated with Feynman 
paths that cover the whole sample. So, we may write 

M 

t ab = J2 A ab(i), (3-17) 

i=l 
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in which A a b(i) is the probability amplitude due to the i-th path connecting the 
channel a to the channel b. We assume that A a i(i) are independent random quater- 
nion variables uniformly distributed on the group U(2). Inserting (3.17) into (2.6) 
we find 

(T a& ) = ^(|A a6 (z)| 2 ), (3.18) 



while 



<^>4£<i^(0iV«60')i s 



4 

^\(T ab )\ (3.19) 

where terms of order unit compared with M were neglected. Using (3.19) we easily 
find 

Var(T ab ) = ^(T a& > 2 , (3.20) 

in agreement with (3.16). 

The second term in (3.13) is a long-range correlation which becomes dominant 
when we consider the fluctuations of the total transmission coefficient = Y2 a Tab 
to the right due to a single mode b excited on the left. From (3.13) we get 

Var ™ = szr.- (321 » 

When a = 1 (conventional waves), (3.21) is the same result obtained in Refs. 5,6 
and 19. In the symplectic case (a = 4), the present model predicts suppression of 
long-range correlations by the same factor 4 found in (3.16). 

The third term in (3.13) describes a uniform long-range correlation between 
all channels. It is the dominant fluctuation of the total transmission coefficient T, 
which is related to the conductance by the Landauer formula [19] G = 2T. From 
(3.13) we find 

Var(T) = I l f , (3.22) 
in agreement with microscopic calculations [11-15]. 
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This result shows that the long-range nature of the diffusive process induces 
non-trivial correlations between transmission coefficients, which in turn generates 
universal fluctuations as observed in the conductance of phase coherent conductors 
[7-10]. 

Let us now consider (3.14). The first term gives for (3 = 1 and (3 = 2 

Var(lU) = (Rab) 2 , (3.23) 

which is the familiar speckle pattern for the reflection coefficient of conventional 
waves. In the symplectic case ((3 = 4), we get 

Var(lU) = ^(iU> 2 , for a ^ 6, (3.24) 

and 

Vav(R aa ) = (R aa ) 2 , (3.25) 

Equation (3.24) can be easily derived following the same reasoning used to obtain 
(3.20). We shall consider the more subtle result shown in (3.25). 

The Feynman paths that contribute to r aa are closed, which implies that inter- 
ference between time reversed trajectories plays an important role and gives rise to 
coherent backscattering. We include this information in r aa by writting it as 

M 

r aa = J2 A aa(*)A aa (i), (3.26) 
i=i 

in which A aa (i) is the probability amplitude associated with the i-th Feynman path 
that starts and finishes at channel a, whilst A aa (i) is the probability amplitude due to 
the time- reversed path. We assume again that A aa (i) are random and independent 
quaternion variables and that the number of paths M is large. Inserting (3.26) into 
(2.5) we find 

<iU> = £(|A aa (z)| 4 >, (3.27) 

i 

while 

(0 = 2^(|A aa (,)| 4 |A aa ( J )| 4 ) 

ij 

^ 2(R aa ) 2 , (3.28) 
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which of course implies (3.25). 

It is very instructive to observe that spin-orbit effects on local fluctuations, 
whilst being intrinsically of quantum mechanical nature, can be understood following 
the same kind of reasoning that explains the familiar speckle pattern of conventional 

waves. 

Another noteworthy feature of (3.14) is that it shows long-range correlations 
which are of the same order of magnitude as those in (3.13). This fact was discussed 
in Ref. 20 to contrast with earlier hypothesis [2], which suggested that absence 
of long-range correlations between reflection coefficients might be consistent with 
universal conductance fluctuations. One can see from (3.14) that it is precisely 
these long-range correlations that give rise to universal fluctuations. We find 



Var(i?) = 




(3.29) 



which is the same value found in (3.22). This is consistent with the identity R = 
N -T. 

Consider now Eq. (3.15). Observe that it shows only long-range correlations, 
since reflected and transmitted waves emerge at different sides of the sample. Such 
correlations exist simply because the interference pattern generated by the scattered 
wave is highly sensitive to any small change in the impurity configuration [23]. 

When a single channel b is excited on the left, we find that the total transmission 
coefficient to the right T b is correlated with the total reflection coefficient to the left 
R b with covariance given by 

COT < T ^> = -3^V? (3 ' 30) 
Comparing (3.30) with (3.21) we find 

Cov(T b ,R b ) = -Var(T 6 ), (3.31) 

which is consistent with the identity R b + T b = 1. 

Finally, the second term in (3.15) is a uniform background correlation that gives 
rise to a universal correlation between the total transmission coefficient T and the 
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total reflection coefficient R. Using (3.12)and (3.15) we find 



Cov(i2, T) = 




(3.32) 



which is consistent with (3.22) and the identity R + T = N. 

IV. Summary and Conclusions 

In this work we have studied the correlations between transmission and reflec- 
tion coefficients of coherent electronic waves diffusively scattered in a disordered 
quasi-one-dimensional medium both with and without spin-orbit scattering. In the 
absence of spin-orbit scattering the quantum mechanical nature of the phenomena 
enters only via the wave-like behavior of the electron and thus the same statisti- 
cal properties of transmitted and reflected intensities are expected to occur in any 
system in which waves propagate coherently through a disordered medium. In the 
presence of spin-orbit scattering, however, new purely quantum mechanical effects 
take place and special statistical properties arise. 

We have obtained novel exact expressions for several correlation functions be- 
tween transmission and reflection coefficients in terms of the number N of eigen- 
channels and certain averages to be evaluated by solving a given evolution equation, 
in which the length L of the sample plays the role of the diffusion time. 

In the metallic regime, which is defined as the one in which the length of the 
sample is much larger than the mean free path, but much smaller than the localiza- 
tion length, we have found explicit expressions for the various correlation functions 
between the reflection and transmission coefficients valid for the unitary, orthogonal 
and symplectic ensembles. We have shown that in the symplectic ensemble spe- 
cial quantum interference effects due to spin-orbit scattering suppresses both local 
fluctuations and long-range correlations of transmission and reflection coefficients. 

The high sensitivity of the interference pattern created by the scattered wave to 
any small change in the impurity configuration leads to interesting correlations be- 
tween the total reflection coefficient and the total transmission coefficient, evaluated 
at opposite sides of the sample, when a single wave mode is excited. We believe that 
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this effect is unique to the diffusive regime and plays an important role in providing 
an understanding of universal fluctuations of global quantities. 

We conclude by remarking that the metallic or diffusive regime is very special in 
that the interferences of the wave fields of the diffusion modes gives rise to remark- 
able statistical properties of transmission and reflection coefficients some of which 
are peculiar to the symplectic ensemble. All long-range correlations that we have 
found in crossed second moments of transmission and reflection coefficients have the 
same order of magnitude and are ultimately responsable for universal fluctuations 
of global quantities like total reflection and total transmission coefficients. The lat- 
ter corresponds to the familiar universal conductance fluctuations observed in phase 
coherent metallic samples. 
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